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Abstract 

In this paper we introduce a new mathematical tool to solve frac- 
tional equations representing models of fractional systems : The Ul- 
* This work was partially supported by Consejo Nacional de Investigaciones Cientificas 

Argentina. 
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tr adistr ibutions . 

Ultradistributions permit us to unify the notion of integral and deriva- 
tive in one only operation. Several examples of application of the 
results obtained are given. 

PACS: 03.65.-W, 03.65.Bz, 03.65.Ca, 03.65.Db. 



2 



1 Introduction 

The use of fractional calculus for modelling physical systems has been consid- 
ered in many works. See for example [HEllS]. We can find also works dealing 
with the application of this mathematical tool in control theory [H [5], [6], [7] . . 

Moreover, there are many physical systems that can be described by 
means of a fractional calculus. Some examples are: chaos [S], long electric 
lines electrochemical process [lO] and dielectric polarization [TT] . 

In this paper we want to introduce a new mathematical framework to 
solve fractional equations representing models of fractional systems which 
was not treated in none of the previous works: The Ultradistributions. 

The paper is organized as follow: in section 2 we introduce definition of 
fractional derivation and integration. In section 3 we give some examples of 
application of the formulae of section 2 using the Fourier Transform and the 
one-side Laplace Transform. In section 3 we present a circuital application. 
Finally in section 4 we discuss the results obtained in sections 1,2 and 3. 
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2 Fractional Calculus 



The purpose of this sections is to introduce definition of fractional derivation 



and integration given in ref . [12] • This definition unifies the notion of integral 
and derivative in one only operation. Let f (x) a distribution of exponential 
type and F(Q) the complex Fourier transformed Tempered Ultradistribution. 
Then: 



¥[C>)=U[3[D.)] 



oo 



f(x)e''''' dx-U[-3{CL]] 



r(x)e'^'' dx (2.1) 



— oo 



(U(x) is the Heaviside step function) and 



f(x) = -L(l)F(a)e-'"'' da 
In 



(2.2) 



where the contour V surround all singularities of F(I1) and runs parallel to 
real axis from — oo to oo above the real axis and from oo to — oo below the 



real axis. According to [12] the fractional derivative of f (x) is given by 



d^f(x) 



1 



dx^ In 



o[-]Q.]^¥[Q.)e-'^'' da + ()(-ja]^a(a)e-'^'' da (2.3) 



Where a(a) is entire analytic and rapidly decreasing. If A = — 1, d^/dx^ is 
the inverse of the derivative (an integration). In this case the second term of 
the right side of (12. 3p gives a primitive of f(x). Using Cauchy's theorem the 
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additional term is 



rafO] 



e-'^^dO = 27ra(0) 



(2.4) 



Of course, an integration should give a primitive plus an arbitrary constant. 
Analogously when A = —2 (a double iterated integration) we have 

aid] 



()■ 



-'^"dO = y + 6x 



(2.5) 



where y and 6 are arbitrary constants. With the change of variables s = — jQ 
formulae (12. ip and (12.21) can be writen as: 



G(s) =U[^H(s) 



f(x)e-^'' dx-U[-9^(s) 



flxje-^" dx 



(2.6) 



and 



Ini 



d) G(s)e='' ds 
r 



(2.7) 



where the contour T surround all singularities of G(S) and runs parallel to 
imaginary axis from — jc» to joo to the right of the imaginary axis and from 
joo to — joo to the left of the imaginary axis. Formula (12. 6p represents the 
two-sided Lapnace Transform. The fractional derivative is now: 



d^ffxl 1 

^Vr^ = — (l)s^G(s)e^'' ds + cbs^alsje^" ds 



dx'*' 2ni ^ 
r 

For the one-side Laplace Transform we have 



(2.8) 



G(s) =U[lH(s)] 



f(x]e-''' dx 



(2.9) 



G(s)e*'' ds 



(2.10) 



a— 500 



and for the fractional derivative: 



d^f(x) _ _1_ 
dx-** In] 



Q+joO 



s^G(s]e''' ds 



(2.11) 



a— joo 



3 Examples 

In this section we give some examples of the application of formulae of the 
precedent section. At first using the Fourier Transform and at second place 
using the one-side Laplace Transform. 

The Fourier Transform 



Let U(x) be the Heaviside step function. 

f(x]=U(x) ; F(a]=up(a)] 

The fractional derivative is: 



e tlx- ^ 



(3.1) 



d^U(x) je-^ 



dx^ 



In 



oupiajja^-^e-^"" da + oD^a[a)e-'^'' dO 



— ]n\ 

In 



(o) + jO^-^e^"'" dcu + oQ^a[D.)e-'^'' dO (3.2) 

r 
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With the use of the result (see ref . [T3] ) 

[w + jO)^-^-^'"'' da) = -In] 



IttA 



r(i-A) 



we obtain: 



d^U(x) 

dx^ ~ r(l -A] 



+ ()a^a(a)e-'"'' dO 



When A = n 



X 



-A 



r(i-A) 



5(n-l) 



()a"a(a]e-'^''da = 



and we have the ordinary derivative: 

d'^U(x) 



dx^ 



When A = —n 



d-^U(x) x'; 



dx 



n! 



+ Qo + aix + aix + ■ • ■ + an-ix 



which is a n-times iterated integral. 

Let 6(x) the Dirac's delta distribution. For it we have: 

f(x) = 5(x) ; F(a) = 

The fractional derivative is: 



d^5(x) _ 
dx^ ~ r(-A) 



When A = n: 



(3.11) 



and when A — — n: 



+ +ao + aix+a2X^ + --- + an-ix^"^ (3.12) 



dx 



[n-1]! 



Let us consider now the fractional derivative of e'^" 



a + b 



(3.13) 



We have: 



dx^ 27T 



a + b 



da + (j)n^a(a)e-'"'' da 



(3.14) 



le 2 
In 



cu + b + jO 



le 2 
In 



cu + b - jO 



(ba^aiaje-'^'' da 



(3.15) 



From ref. [T^ we obtain: 



(x + y]^ _ 



x+ |3 



27rU(p) ^^^ ,, P'VP^ct)[-A,1-A,j(y-|3)p] 



(3.16) 
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where cj) is the confluent hypergeometric function. Thus the fractional deriva- 



tive is: 



dx^ r(i-A) 

With the use of equality: 



cl)(1,1 -AJbx) +()a^a(a]e"'^'' dQ (3.17) 



(t)(1 , 1 - A, jbx) = (jbx) V^'' [r(1 - A] + Ar(-A, jbx) 



(3.18) 



where r(zi,Z2) is the incomplete gamma function, f l3.17p takes the form: 



dx^ 



A jbx 



1 + -F7^ — -r(-Ajbx) 



r(i-A) 



+ 



When X — n 



dx^ 



n^jbx 



and when A — — n: 

^-rigjbx 



dx" 



TL-1 



(3.19) 



(3.20) 



(3.21) 



The Laplace Transform 



If we use the one-side Laplace transform to evaluate the fractional derivative 



of U(x),then: 

f(x)=U(x) ; G(s)=U[9^(s) 



(3.22) 
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and as a consequence: 



d^U(x) _ _1_ 



a+)oo 



U[9^(s]]s^-^e^'' ds 



a—) 00 



In 



.-A 



[a + ]s 



l-A 



ds 



r(i-A) 



.-A 



d^U(x) 

dx^ ~ r(1 -A] 



When A = n we obtain 



d^U(x) 
dx^ 



5(n-l) 



which coincides with f l3.7l) . When A = — n the resuh is: 



d-^U(x) _ x^ 
dx-^ ^ TlJ 



In a analog way we obtain for Dirac's delta distribution: 



d^6(x) 




dx^ 


r(-A] 


d^5(x) 




dx^ 




d-^6(x) 


^+ 


dx-^ 


(n-ll! 



Finally we consuder the fractional derivative of e' : 

f(x)=U(x)e'^'' ; G(s) = HM^ 

s — 1 b 
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According to fl2.11l) : 



dx^ In] 



a+)oo 



s — jb 



a— )oo 



e 2 
2n] 



s + b + jO^ 



And thus: 



d^Ufxle'^'' Ufxlx-^ 



dx^ r(l-A) 
Using fl3.18p . (13.341) transforms into: 



-AJbx) 



d^Ulxje'^" 
dx^ 



= (jb)^U(x)e 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



When A = n: 



and when A — — n: 



dx^ 



dx-^ 



(jb)^U(x)e'^'' 



(3.36) 



(3.37) 



4 Circuital Application 

As circuital apphcation we consider a semi-infinite cable with a voltage V = 
Voc''"* applied at one end. We use first the Fourier transform and then the 
Laplace transform for see the diferences between both treatments. 
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The Fourier Transform 



We should solve the system: 



a2f(x,t) ^^af(x,t) 



at 



= ; X > 



(4.1) 



f(0,t) =Voe'-* 

where R is the resistance per unit length and C is the capacitance per unit 
length. Let V(x, t) the voltage along the semi-infinite cable. We use a for- 
malism developed in ref. [15j to solve the system (14. ip . It consist in to define: 



V(x,t) =U(x)f(x,t) 

af(x,t) 



ax 



x=0 



The differential equation in fl4.ip transforms into: 



RC — — — = (xjVoe' + 6(xjg(t] 



at 



(4.2) 



(4.3) 



Taking the Fourier transform of (14.30 we obtain: 



'^i(x-\, ocz) = 7rjVo6(a.i -|- cv] 



i-i. 



OCz + ^yZ—OC^RC 



(1 -j)V-2cciRC 



(4.4) 
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1 1 

Oil — ^-^\/—OL\^Q- 0L2 + ^A/-aiRC 
Deprecating the exponential increasing in the solution we obtain: 

^(ai) = -(1 + j)7rv/-2cxiRC b[oi^ + w] 



and then we obtain: 



g(t) 



The current i(x, t) is: 



i(x, t) 



1 3V(x,t) 

R ax 



; X > 



As: 



3V(x,t) 
dx 



(1+3) 



<^RC, , / oi R C „ 



then: 



i(x,t) = (1 + jj^/^Voe-v^V'-^-v^"' ; X > 



If we take A = 1 /2 in ( 13.191 we obtain: 



2V7r 2 



+ ()Z^a(Z)e-^^*dZ 
r 
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Ot2 



1 r., 1 . 

2V7r 2 



(|)Zia(Z,x)e-'^MZ (4.13) 



Thus we have a relation between the current and the time derivative of the 



voltage: 



i(x, t] 



(jcu)7r(-ijcut) 



V(x,t] 



()Zia(Z,x)e-'^MZ 



(4.14) 



If we consider only the first term in the rigth side of fl4.14p we obtain the 
more habitual result: 

i(x,t| = A/— , (4.15) 



R dti 



The Laplace Transform 

If we use the Laplace transform in place of the Fourier transform to evaluate 
the fractional derivatives, (14. 121) . (14. 131) and (14.141) are replaced by: 



= (ja))ie'-* 

dt2 



1 + 7^r(-l jcut) 

2^/7r 2 



(4.16) 



= 



3t2 



2^/7l 2 



■r ^ / C 

i(x,tj = \/-^ 



3i (ja;)^r(-^,ja)t) 



dti 



2v^ 



V(x,t) 



(4.18) 



14 



Difference between tfiis results and tfie precedents is tfie term tliat contain 
a contour integral. 

5 Discussion 

In this paper we have shown that Ultradistribution Theory is an adequate 
framework to define a Fractional Caculus and its applications. This definition 
unifies the notion of integral and derivative in one only operation. Several 
examples of application of fractional derivative are given, including a circuital 
application: a semi-infinite cable with a voltage V = Voe^"'* applied at one 
end. 
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